Prom local to nonlocal Fermi liquid in doped antiferromagnets 
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The variation of single-particle spectral functions with doping is studied numerically within the 
t-J model. It is shown that corresponding self energies change from local ones at the intermediate 
doping to strongly nonlocal ones for a weakly doped antiferromagnet. The nonlocality shows up 
most clearly in the pseudogap emerging in the density of states, due to the onset of short-range 
antiferromagnetic correlations. 
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The understanding of electron spectral properties in 
the normal-state of cuprate superconductors remains one 
of the hardest challenges for theoreticians. The single- 
electron spectral function is directly related to angle- 
resolved photoemission spectroscopy (ARPES) experi- 
ments, which made impressive progress in the last decade 
jij, in particular in the studies of the development as 
function of the hole-doping [|||] with respect to the ref- 
erence antiferromagnetic (AFM) material. ARPES on 
optimally doped materials reveals electronic excitations 
with a large Fermi surface (FS) with the volume con- 
sistent with the normal Fermi liquid (FL) picture. The 
shape of spectral functions is more controversial [1 , being 
possibly closer to the concept of a marginal FL [1 than 
that of a usual Landau FL. Recently most interesting 
ARPES results have been obtained for underdoped Bi- 
SrCaCuO compounds [^,|| , where a leading-edge shift at 
T > T c indicates a feature consistent with d-wave super- 
conductivity persisting in the normal state. A pseudogap 
(although much larger) is observed also in the integrated 
photoemission of LaSrCuO cuprates || . These phenom- 
ena are closely related to relevant energy scales in doped 
AFM. In underdoped cuprates experiments confirm that 
besides T c a crossover T* > T c |J, pronounced in the 
spin susceptibility xo an d in transport quantities p(T), 
Rh(T), and a lower one T sg , most clearly seen in the 
NMR relaxation. 

Most theoretical calculations of spectral functions 
A(k, to) in doped 2D AFM rely on numerical approaches 
[7|. The exact diagonalization results for the t-J model 
[3 1 and the quantum Monte Carlo results for the Hub- 
bard model H provide the evidence that prototype mod- 
els indeed yield a large FS at the intermediate doping 
with the quasiparticle (QP) dispersion close to the exper- 
iment and a 2D tight-binding band. A recent application 
of the finite-temperature Lanczos method jl(],[yj in the 
same regime allows for additional conclusions for spectral 
functions Jl^ ]: a) A(\t, us) is strongly asymmetric with re- 
spect to the FS, with underdamped QP for k > kp and 
essentially overdamped QP for k < kp, c) the self-energy 



follows the marginal FL form ImE ~ j(\ui\ + £T). The 
regime of low but finite doping is even harder to exam- 
ine. Among the few studies of A(k, uS) only the quantum 
Monte Carlo results jl3| seem to indicate the existence 
of a pseudogap. 

In the present Letter we point out two new and unex- 
pected features of the t-J model as a function of doping. 
The first is that at intermediate doping A(k, u>) is well ac- 
counted for by the assumption of the local FL |lj] , where 
the self-energy is like that recently invoked in the context 
of infinite-dimensional models p| . The second new fea- 
ture is in the underdoped regime. Reducing the doping 
towards an undoped AFM the nonlocality becomes essen- 
tial, leading to a pseudogap in the single-electron density 
of states (DOS). 

We study the planar t-J model as the prototype model 
for strongly correlated electrons in a doped AFM, 
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with projected operators Ci S = Ci S (l — rij,- s ) preventing 
the double occupation of sites. The Green's function can 
be defined in the usual way G(k, u>) — ((5k s ; cj^})^^ 
where fj, is the chemical potential. 

It should be stressed that A(k, u) = — (l/-7r)ImG(k, w) 
is not normalized to unity |8[|l^], since ({ck s , 5^}+) = 
(1 + c/,)/2 = a < 1. This introduces some ambiguity 
in the definition of S. We choose here the one which 
guarantees the analytical behavior, S(k, u) — ► ±oo) = 
p6[ (in contrast to the definition in Ref. Jl2[), 
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The 'free' propagation term (k is nontrivial, since the 
hopping term in Eq.(Q) involves projected operators. It 
is given by [[ill 
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where 7k = J2j nn o ex P(*k ' r j)- Note that 77 represents 
the overall renormalization of the spectral dispersion, not 
just of the QP part. Since n.n. correlation (SfSJ) de- 
creases with doping one can conclude from Eq.(0) that 
at the intermediate doping 77 ~ a < 0.5 while in a weakly 
doped AFM ri ~ 0.27. 

We focus on the doping dependence of the A(k, uj), the 
corresponding £(k, uj) and the DOS, defined by j\f(e) = 
(2/iV) ^ k A(k, e — fi). Calculations are performed at 
T > and on small square lattices with N = 16, 18 and 
20 sites using the finite- T Lanczos method [10 llj]. We 
fix J = 0.3 t to be representative of strongly correlated 
systems and of cuprates in particular, where t ~ 0.4 eV. 
It is important to realize that the results are meaningful 
for T > Tt s . Tf s is the temperature where finite-size ef- 
fects start to dominate. In available systems Tf s > 0.1 t 
at the intermediate doping 0.1 < Ch < 0.25 while it is 
increasing towards undoped AFM, e.g. Tj s > 0.2 t for 

Results for A(k, uj) at intermediate doping have been 
partly presented and discussed previously [jl2j . We point 
out here a novel aspect of the results, i.e. a nearly k 
independent self-energy, S(k, uj) ~ El(oj). We show in 
Fig. la ReS obtained for Ch = 3/16, which reveal a weak 
variation with k. This behavior is characterstic of a local 
FL |TJ,|l^]. This is in a sharp contrast to a strong k 
dependence of the self-energy in the underdoped case in 
Fig. lb. One of the essential consequences of the local FL 
is that the DOS at the FS A/"(m) is independent of the self- 
energy, being (usually) equal to that of noninteracting 
fcrmions. Note, however, that in the t-J model also the 
'free' electron dispersion (k is renormalized, as follows 
from Eq.(|). Nevertheless the DOS for c h = 3/16,4/16 
p2| reveal a smooth variation j\f(e ~ (m) oc 1 /a. 
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The local FL character is lost entering the underdoped 
regime. We show in Fig. 2 A(k, uj) for the lowest concen- 
trations Ch — and ~ 0.06, combining systems with 
N = 16 and N = 18 sites, for the lowest T ~T fs . Rather 
well understood limit (at least for T — > 0) is Ch = 0, 
where one is studying a hole in an spin background with 
the long range AFM order. Note that T > and Ch — > 
requires /j, — > —00 hence it is meaningful to present 
A(k, e). Consistent with analytical approximations JTt| ] 
the spectral function shows a coherent QP peak with a 
dispersion on the energy scale determined by J (band- 
width W ~ 2 J), with the maximum at k* = (n/2,n/2) 
with eo ~ 2t. In addition to a QP feature there is a 
broad incoherent background for e <C eo- The latter has 
still some structure, which must be related to the AFM 
order, since it disappears for T > J. For Ch = the 
dispersion corresponds to a double AFM unit cell hence 
the local FL cannot apply. 
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FIG. 1. Self energies Re£(k, u>) for the doping: a) 
c h = 3/16, and b) c h = 1/18, both for T = 0.15 t. 



FIG. 2. Spectral functions for the undoped AFM (left), 
and for system with finite doping ch = l/N, where results for 
N — 16, 18 are combined for fixed T = 0.15 t. 

More challenging are results for c/, ~ 0.06, also shown 
in Fig. 2. A(k, uj) in this regime merge some properties 
of the undoped AFM and of the intermediate doping. As 
expected for > there is besides a hole (uj < 0) part 
also an electron part with uj > 0, its integrated intensity 
increasing as 2c^. It is evident that the electron part 
is significant only for k > kp, i.e. for k outside the FS 
which seems to be close to the square AFM Brillouin 
zone. Another fact is that in A(k > kp,ui) with a QP 
peak at uj > there is also a shadow feature for uj < 0, 
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most pronounced for k = (71", 7r/3), (ir, n/2) as well as for 
k = Q = (tt, 7r). In fact the hole part shows a dispersion 
very analogous to the one at Ch — 0, since the dispersion 
seems to fold back for k > kp. This is best seen near the 
X point k = (ir,0) where the hole dispersion does not 
seem to reach uj = 0, hence an apparent gap between the 
hole and the electron QP persists. The folding effect is 
hardly visible along the diagonal T-M, i.e. for k = (x, x), 
where the QP peak disperses through the FS in a more 
normal way. 

Gap features appear best visible in the DOS which 
as a local quantity is less affected by finite size effects. 
In Fig. 3 we show j\f(e) for three lowest dopings Ch = 
0/20,1/20,2/18 and different T < J. For c h > the 
vertical line denotes the chemical potential fi(T = 0.1 t). 
While at larger doping the DOS does not show any sign 
of a gap at the the Fermi energy e ~ /1 (at least not for 
T > Tf s ) it is evident from Fig. 3 that the pseudogap 
starts to emerge for Ch = 2/18 and is well pronounced for 
c h = 1/20. 
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FIG. 3. Density of states J\f(e) for systems with various 
doping c h = 0/18, 1/18, 2/18 and different T. The thin verti- 
cal lines denote the chemical potemtial fi(T = 0.1 t). 

The temperature dependence of the DOS is very in- 
structive. It is plausible that for Ch — the QP (spin- 
polaron) propagation (Tt]] becomes incoherent for T > J 
due to lost AFM short range order, hence a QP peak in 
j\f(e) broadens. It is interesting to note that at T ~ J 
some sharper resonances as well disappear in the 'inco- 
herent' part e -C eo- These resonant features at Ch. = 
and T — > can be interpreted as real excited resonance 



states of the spin-polaron |l7|j . The latter are smeared 
out due to the spin disorder at T > J or due to the pres- 
ence of other holes as concluded from Fig. 3 for c h > 
where no such resonances are visible. 

A novel feature is the pseudogap at e ~ p,. The latter 
is not a finite-size effect since very similar results are ob- 
tained e.g. for systems 1/16 and 1/18 JE^ |. The structure 
disappears for T > J, hence it must be related to the on- 
set of the short-range AFM order. Also its energy span 
Ae < t seems to be determined rather by J. The pseudo- 
gap is inconsistent with the concept of the local FL, since 
the latter assumption would lead to an enhanced jV(/j.) 
for Ch — > due to reduced 77, Eq.(||). Hence the pseudo- 
gap can be accounted only by a pronounced nonlocality, 
i.e. k dependence of E(k ~ kp , cj ~ 0). 

The underlying scenario is that the system remains a 
FL with a well defined FS and a conserved (large) FS vol- 
ume. This requires at least ImE(ki?,0) — > for T — > 
whereby the FS is determined by (k F = A* ~ ReE(kp, 0). 
Our numerical results for a weakly doped AFM seem to 
be qualitatively consistent with these assumptions. Note 
that alternative scenarios have been proposed for the 
regime of weak doping e.g. invoking a small hole-pocket 
FS The FL assumption leads to the usual form for 
the DOS, 
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The pseudogap thus requires a large \dRe E(k, 0)/<ik at 
k = kp. In Fig. lb we show Re S(k, uj) at Ch — 1/18 and 
fixed T = 0.15 t. First we realize that ReS(k, 0) = /i — Ck 
locates the presumed FS again close to the AFM zone 
boundary. The differences with respect to intermediate 
doping are easily observed. While for k inside the 
FS the self-energy is quite similar to the local FL, there 
is an increasing deviation from the latter for k > kp. In 
particular this shows up as an oscillation in Re£(k,a->) 
most pronounced for k = Q. The effect is so strong 
that it leads for most fc > fcp to a double solution for a 
QP peak, i.e. i?k = Ck + ReS(k, Et) — fx, explaining the 
shadow features and folded spectra in Fig. 2 for most k 
outside the FS. 

The double solutions for are analogous to those 
appearing within the spin-bag theory |l9[ taking into ac- 
count the AFM short-range order. The essential differ- 
ence is however an evident asymmetry since no such effect 
appears for k inside the FS. To simulate the observed cor- 
rections to El due to AFM correlations can be designed. 
For a long-range AFM order one could expect the form 



E(k, w) - £ L H + o(k)G(k + Q, uj). 



(5) 



The expression (^) naturally leads to oscillations in 
E(k, uj) and possible double QP solutions, in particular 
near the half-filling where the FS is close to the nest- 
ing AFM zone boundary. In order to account for the 
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pronounced asymmetry with respect to FS one has to as- 
sume an essential dependence of a(k), i.e. a(k <C fcf) ~ 
and its strong increase near ftp. 

It is however still feasible that in a weakly doped AFM 
the emerging pseudogap coexists with a well defined large 
FS. AFM correlations open a pseudogap between the 
hole-like and electron-like QP resonances coexisting for 
k ~ \ip , visible in particular near the X point in Fig. 2. 
Here it seems that the QP dispersion does not traverse 
the FS. But the assumption ImS(k, 0) = at T -> 
leads to another QP crossing the FS. The latter peak 
must be of a small QP weight Z = Z^ F . This could be 
consistent with the inequality Z < 2ch [fl2f derived on 
the assumption of the monotonous fall-off of n^. Such a 
QP excitation would be hard to observe in numerical cal- 
culations as well as in experiments since it can be easily 
overdamped either by T > or other effects. 

Let us turn to the discussion of the relation of our 
results to photoemission experiments on cuprates. The 
latter capture mostly the hole-part lu < 0. For the in- 
termediate doping our A(k, uS) are qualitatively consis- 
tent with ARPES measurements , as discussed in Ref. 
[13. ARPES on underdoped BiSrCaCuO compounds 
\W shows an opening of a pseudogap near the X point 
where the FS seems to disappear. On the other hand 
the QP peak disperses through the FS along the T-M di- 
rection, rather like in a normal FL. Both facts are well 
consistent with results in Fig. 2. Recently the DOS has 
been measured via the integrated photoemission (also the 
inverse part) for LaSrCuO material in a wide range of 
doping [El. More reliable are emission spectra related to 
N + {lu) — J\[(u> + n)f(uj), giving information for lu < 0. 
In the overdoped systems Ch > 0.2 N + {lu < 0) appears 
flat. On the other hand, for Ch < 0.17 a pseudogap starts 
to emerge gradually, e.g. the inflection point in A/" + (w) 
moves from lu ~ towards u> ~ — 0.2eV at lowest doping. 
This is quite close to our results and values in Fig. 3. The 
pseudogap scale and its doping dependence seems thus to 
be related to T* ||, appearing e.g. as the maximum in 
Xo- Note however that such a gap is much larger than 
the leading-edge shift - 30meV found in ARPES @ in 
underdoped samples. 

In conclusion, we have shown that within the t-J model 
the spectral function changes with doping from that of a 
local FL to a strongly nonlocal one. It seems plausible 
that the locality of the self-energy at intermediate doping 
is related to the coupling of QP to large (mostly local) 
quantum fluctuations of the spin subsystem [jl6| whereby 
its anomalous behavior emerges due to the frustration in- 
duced by holes. The locality has several important conse- 
quences. Apart form the mentioned effect on j\f(fJ,), one 
key feature is the reduction of FL parameters to only two, 
which are even simply related It is still a question 
whether such a local FL has properties analogous to ones 
obtained in infinite-D models [^5) . 

In contrast, the nonlocality of £(k,w) induced by 



short-range AFM correlations becomes essential for the 
regime of weak doping, and is dominating the low-energy 
behavior lu ~ 0. The most pronounced effect is the open- 
ing of the pseudogap in the DOS. A(k, lu) in this regime 
can reveal a coexistence of a hole-like dispersion anal- 
ogous to undoped AFM and an electron-like resonance. 
One of the challenging questions remains whether in un- 
derdoped systems we are still dealing with a well defined 
large FS, whereby our results do not contradict this sce- 
nario. It should also be pointed out that nonlocality is 
crucial to account for the marked difference in the low- 
doping regime between nearly free-fermion-like Wilson 
ratio and a strongly enhanced (possibly even singular) 
compressibility fllJl , indicating the closeness to a charge 
instability. 
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